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ABSTRACT 

Using hydrodynamic simulations and a technique to extract the rotational component of the velocity field, 
we show how bubbles of relativistic gas inflated by AGN jets in galaxy clusters act as a catalyst, transforming 
the energy carried by sound and shock waves to heat. The energy is stored in a vortex field around the bubbles 
which can subsequently be dissipated. The efficiency of this process is set mainly by the fraction of the cluster 
volume filled by (sub-)kpc scale filaments and bubbles of relativistic plasma. 

Subject headings: hydrodynamics — instabilities — shock waves — methods: numerical — galaxies: clusters: 
general — ISM: bubbles 



1. INTRODUCTION 

Models for the formation of galaxy clusters require an ex- 
ternal source of energy to avoid catastrophic cooling of gas at 
the centers of massive clusters and to explain the observed 
lack of cool gas ( Peterson et al. 2003) in so-called cooling 
flow clusters (.Fabian_199 4). Relativistic jets from supermas- 
sive black holes are th e most likely source of this energy 
iBinnev & Taboiiri995ft . While jets may carry sufficient ki- 
netic energy, it has been unclear how efficiently that energy 
can be transferred to heat the cluster gas. 

When jets interact with the gas in galaxy clusters (the intra- 
cluster medium, ICM), they invariably inflate hot, tenuous 
bubbles of relativistic gas, observed as synchrotron radio neb- 
ulae (Scheuer 1974), a.k.a. radio lobes. Evidence is mounting 
(|En6Hn 1999; Birzan et al. 2004; Giovannini & Feretti 2004) 
that the ICM is littered with such bubbles, each of which 
starts out growi ng supersonically, driving a strong shock 
JHeinz et alJfT99 8) that efficiently heats the ICM it encoun- 
ters. However, X-ray observations show that the expansion 
rapidly decelerates and most bubbles expa nd subson ically, 
therefore not producing strong shocks (tFabian et alJ i200Q: 
iMcNamara et al. 2000). 

The bubbles themselves store a fraction of the initial jet en- 
ergy, some of which can be transferred to the gas and dis- 
sipated duri ng the bubble's slow, buoyant rise in t he cluster 
atmosphere JChurazov et alJ200lHBegel mani200l1). Another 
significant fraction of the jet en ergy is released in the form 
of sound and weak shock waves ("Fabia n et alii200 3a). In ho- 
mogeneous media, dissipation of wave energy is determined 
by microscopic transport processes like viscosity or conduc- 
tion, which are well known for unmagnetized plasmas. Re- 
alistically, though, the presence of even weak i3-fields in the 
ICM can alter the transport coefficients, making estimates of 
the dissipation rate very uncertain. High values for the micro- 
scopic viscosity have been post ulated as a way to d issipate the 
wave energy t,Ruszkowski et aL.2004; iRevnolds et aL.2005c 
iFabian et alJ2005ft . 

In this letter we argue that the presence of filaments and 
bubbles of relativistic gas provides an efficient way to extract 
the wave energy and heat the ICM even //viscosity and ther- 
mal conduction are strongly suppressed. In fjJlwe review the 
process responsible for this dissipation, ^ presents the nu- 



FlG. 1 . — Map of fluid density (color scale) and velocity field (arrows) for 
Mach 1 .9 shock running over a bubble with density contrast of 1 : 1 00. Shown 
are time frames and 31 in units of the sound crossing time of the bubble. 



merical methods we used, [j^discusses how our results can be 
applied to cluster heating, and i|5lsummarizes. 

2. THE RICHTMYER-MESHKOV INSTABILITY 

The energy transfer to the ICM occurs through a pro- 
cess known as th e Richtmy er-Mes hkov instability (RMI 
'Richtmver 196^; 'Meshkovl fT969t lOuirk & Karnil [l995 
Inogamov 1999). A plane shock or sound wave passing over a 
density discontinuity that is inclined with respect to the wave 
front, as it occurs at the boundary of the ICM with a rela- 
tivistic bub ble, induces a well-loca lized vortex flow around 
the bubble (EnBlin & Briiggen 2002), as illustrated in Fig. 1. 
A qualitatively similar process occurs when a shock passes 
over o verdense clouds iKlein et anil994HKornreich & Scald 
120001) . 

Because the high sound speed keeps the material inside the 
bubble at uniform (but time varying) pressure, the bubble acts 
like a hydraulic piston, distributing pressure applied at one 
end over its entire surface. The arrival of the wave sets up 
a pressure gradient which accelerates the bubble/wave inter- 
face inward, introducing significant inward transverse motion 
and distorting the bubble. Ahead of the wave, the overpres- 
sured bubble expands into the ICM, inducing significant out- 
ward transverse motion. Once the shock front has passed over 
the bubble, it leaves behind a vortex ring in what was pre- 
viously irrotational flow. For a compression wave of finite 
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width A, the passage of the traihng end of the wave will act 
in the opposite sense, partly restoring the bubble's original 
shape and canceling previously generated vorticity. This re- 
duces the amount of energy deposited by the wave if A is not 
large compared to the bubble radius R. 

In the limit of low bubble density pbubbio ^ PiCM, and 
for R ^ \, one can easily estimate the amount of the energy 
funneled into this vortex: Material of density picM enters the 
bubble volume Vbubbic with velocity wicm (the bulk veloc- 
ity of the compressed gas). Thus, the energy trapped in the 
toroidal vortex should be 

1 2 
-E-rot — 52^bubblcPlCMWiCM (1) 

where g depends on the bubble's geometry and should be of 
order unity for spherical bubbles. This energy is extracted 
from the wave and available for dissipation in the ICM. 

In order to verify the dimensional analysis that leads to 
eq. O and to calibrate the coefficient g, we undertook a set 
of numerical hydrodynamic simulations of shock-bubble in- 
teractions. We shall present the details of our numerical study 
before discussing the results. 

3. NUMERICAL METHOD 

We performed multiple parameter studies of shock-bubble 
interactions. The setup consisted of a hot bubble of varying 
geometry embedded in a colder, denser medium. Shocks or 
non-linear sound waves of varying finite width A and pulse 
height (i.e., Mach number) were introduced to travel across 
the bubble. 

The simulations were carried out using the publically avail- 
able FLASH code (Frvxell et al. 2000). FLASH is an adaptive 
mesh refinement hydrodynamics code that uses a second order 
accurate piece-wise parabolic solver Most of our simulations 
were carried out on a 2D Cartesian grid with an effective grid 
size of 2048x3072 cells and an effective resolution of 64 cells 
across the bubble. 

Figure 1 shows a typical setup with a density contrast of 
100 and pressure balance across the bubble, the ratio of sound 
speeds is 10 from inside to outside of the bubble. We assumed 
an adiabatic equation of state except for the shock, which is 
handled implicitly by a shock capturing scheme. The ratio 
of specific heats for both hot and cold fluids was 7 — 5/3. 
While this will not affect the results qualitatively, a softer 7 = 
4/3 equation of state would increase the compressibility of 
the bubble and lead to a slightly higher value of g. 

In most of the simulations the shock was taken to be a semi- 
infinite piston with a sharp discontinuity in pressure, density 
and veloci ty, satisfying the Rankin e-Hugoniot shock jump 
conditions. lLandau & Lifs hitz'1987, In the case of finite pulse 
width A (Fig. 2), the initial pulse profile was taken to be a top- 
hat function in all fluid variables. 

Our choice of grid implies a cylindrical bubble geometry. 
To confirm that our conclusions are independent of geometry 
we ran test-cases of spherical bubbles on a 2D axi-symmetric 
grid. A resolution study confirmed that our simulations are 
independent of numerical resolution. The results are insensi- 
tive to the bubble/environment density ratio Pbubbio//OiCM as 
long as pbubbic ^ PiCM- While the enforced symmetry of 
our 2D simulations might slightly alter the amount of vortic- 
ity generated and thus the numerical coefficient g, it has been 
shown that the RMI does operate in 3D and the results of this 
paper will most likely not be affected severely by restricting 
the analysis to 2D. 
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Fig. 2. — Black: Kinetic energy in the rotational velocity field ex- 
tracted from sound/shock wave as a function of the bulk velocity of 
the shocked gas, vicm- Points correspond to Mach numbers M = 
1.02, 1.07, 1.22, 1.87, 4.20. Grey: Energy extraction efficiency g as func- 
tion of pulse width A (relative to maximal value g = 0.97 for A — > 00). 

The rotational component of the velocity field was extracted 
by solving the vorticity equation Vxv = Vx(VxA) for 
the vector potential A (which reduces to Poisson's equation 
for —A in 2D Cartesian coordinates), subject to the bound- 
ary condition of vanishing rotational velocity at infinity, mak- 
ing Wrot invariant under Galilean transformation. We used 
a Fast Fourier Transform method to solve the equation nu- 
merically. We then calculated the kinetic energy i^iot = 
/ dVpCV X A)^/2 contained in the rotational part of the flow. 
The potential velocity field can be found by solving Poisson's 
equation = Vv. 

4. DISCUSSION 
4.1. Numerical calibration of g 

These simulations can be used to test eq. Q and to cali- 
brate the numerical coefficient g. Fig. 2 shows the kinetic en- 
ergy ii^iot in rotational flow as a function of uicm in the limit 
R ^ A. The numerical results are in good agreement with the 
estimate from eq. Q, with a numerical fit of g = 0.97. Fig. 2 
also shows that for A < IQR, g is proportional to (A/i?)^. 
This dependence of g on A/i? is linked to the degree of bub- 
ble deformation produced during shock passage. Because the 
total amount of energy in the wave is proportional to A, the en- 
ergy transfer efficiency is 77 = g{R/X) cx \/R. Thus, small 
bubbles are more efficient at extracting energy than large ones. 

Bubble shape and orientation also affect g: Fig. 3 shows 
that i^iot depends on w and I, the bubble's dimensions paral- 
lel and perpendicular to the wave. The numerical fit for as- 
pect ratios l/w not too far from unity is 5 cx (l/w)'^'^: the 
more surface area the bubble presents to the wave relative to 
its volume, the less efficient the mechanism becomes. The de- 
pendence of g on filament orientation shown in Fig. 3 (grey) 
is due to the effective aspect ratio changing with angle. This 
dependence can be easily understood considering that all ef- 
fects essential for deposition of rotational energy into the ICM 
(vorticity generation and bubble deformation) scale with the 
length I of the filament perpendicular to the front. 

4.2. Application to galaxy clusters 
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Fig. 3. — Extraction efficiency g for kinetic energy in rotational flow as 
a function of aspect ratio I / w (black circles) of the filament and best fit re- 
gression oi g oc {l/w)°-^ (black line). Bubbles exposing less surface area 
to the shock are more efficient at extracting energy per unit bubble volume. 
Also shown: Efficiency <; as a function of angle between filament and shock 
normal for an aspect ratio oi l/w = 3.5 (grey squares) and expected depen- 
dence g oc (Ay / AxY^'^ from measured x- and y-cuts across the filament 
(grey Une). 



These resuhs imply that a distribution of bubbles of hot gas 
whose characteristic size is small compared to their distance 
r from the cluster center can extract energy from a sound or 
shock wave passing over them rather efficiently. The total en- 
ergy captured equals the energy density in the wave times the 
volume filled by bubbles. To demonstrate this. Fig. 4 shows a 
simulation of a random distribution of bubbles (filling about 
10% of the volume) before and after the shock passage. The 
top panels show the creation of a turbulence field around the 
bubbles in the wake of the shock. The bottom panel shows 
the decomposition of the velocity field into rotational and po- 
tential flow, i^iot and the viscous dissipation rate are well lo- 
calized around the bubbles even long after the shock passage. 

The velocity in the differentially rotating vortices drops 
with distance from the bubble center and a number of dis- 
sipative processes (such as the magneto-rotational instability 
(jBalbus & Hawlev 1991), turbulence, and microscopic vis- 
cosity) will dissipate the energy contained in the vortex. The 
dissipated heat and entropy will go directly into the thermal 
gas surrounding the bubble. Since small bubbles rise through 
the cluster gas with spe eds much smaller than the sound speed 
dChurazov et alJ l2002h . the deposition of rotational energy 
and the subsequent heating occur at an essentially fixed lo- 
cation in the cluster compared to the passing wave. 

Whether the condition i? <C A is satisfied in galaxy clusters 
depends on the unknown size distribution of bubbles in the 
ICM. In sub-sonic models of ICM heating (Begelmani l200!l 
IChurazov et alJ2 0'()^ it is typically assumed that a large num- 
ber of small bubbles is released rather than few large bub- 
bles. X-ray and radio images of nearby clusters indicate 
that the distributio n of radio pl asma c an be rather complex 
jQwen etalJl20(Kit ^oung et aj] 11001 iFabian et al. 2003a), 
though resolution limits and the line of sight projection in X- 
ray and radio maps make a direct detection of small filaments 
impossible. The bubbles detected directly in these clusters 
are of similar size to the sound waves seen, implying a low 
energy extraction efficiency g. In addition to these prominent. 




Fig. 4. — Random distribution of low density bubbles (filling fraction 
/ = 10%) before (panel 1) and after (panel 2) passage of an initially square 
pulse of Mach 1.87, which induces a complex turbulence field in the wake of 
the bubbles and distoils the wave front, compai'ed to the same setup but with 
uniform density (panel 3). Panel 4: Kinetic energy of the decomposed ve- 
locity field: rotational flow (red); potential flow (sound/shock waves, green); 
and viscous dissipation rate (blue, in units of specific viscosity). Note that 
the rotational energy is localized around regions of enhanced bubble num- 
ber density and that some of the energy in potential flow is radiated away in 
sound waves. 



large radio bubbles, radio observ ations of clusters a lso show 
more amorphous radio halos (e.g. 'P edlar et al .^l990l). extend- 
ing over the central ~ 100 kpc. Such halos could indicate that 
smaller radio bubbles/filaments are present in the ICM. 

Dynamical (Rayleigh-Taylor and Kelvin-Helmholtz) insta- 
bilities ( Landau & Lifshitz 1987) lead to fragmentation and 
shredding of individual bubbles, producing a spectrum of 
smaller filaments. In the absence of magnetic fields, buoy- 
antly rising bubbles are shredded into smaller structures while 
passing a distance comparable to their own diameter, within 
a time scale r^hrod 2 x lO'^ yrs^(i?/10 kpc)(4 keV/T) 
for a cluster of temperature T. As long as Tghrod is com- 
parable to or smaller than the buoyant rise time Tbuoy ^ 

5 X Tshrod[(^cooi/100kpc)(10kpc/i?)]"^^^ a significant frac- 
tion of a bubble's volume will be shredded within the cluster's 
cooling radius rcooi- A more quantitative evaluation of bub- 
ble shredding, including the effects of the magnetic fields and 
viscosity, which can significantly suppress the shredding ef- 
ficiency, is a topic of ongoing research (Revnolds et al. 200^ 
'Kaiser et al."2005^ and beyond the scope of this paper 

Shredding bubbles to smaller sizes increases the efficiency 
g until they reach a characteristic size much smaller than the 
wave length of the perturbations. After this point, g is unaf- 
fected by bubble size, though the buoyant rise velocity still 
decreases with bubble size, thus leaving more and more time 
for dissipation processes to act on the vorticity field. Ulti- 
mately, if bubbles get shredded to small enough sizes their 
relativistic gas content might get microscopically mixed with 
the ICM. The energy contained in the low energy electrons 
can then also contribute to ICM heating via Coulomb losses. 

As a possible application of the process suggested above, 
consider the Perseus cluster, one of the brightest cooling flow 
clusters. The size of the prominent central radio bubbles 
is about 10 kpc. This sets the typical thickness of A ^10 
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kpc of the waves generated by bubble inflatio n, consistent 
with th e concentric ripples seen in X-ray images JFabian et alJ 
120031 . Since the bubbles are comparable in size to A, they 
are inefficient at extracting energy from AGN induced sound 
waves. However, if the ICM in Perseus is filled with numer- 
ous smaller bubbles, the typical attenuation length for a com- 
pression wave will he L ^ A//, where / is the fraction of 
cluster volume filled by bubbles. For / ~ 10%, most of the 
wave energy would be trapped in the inner ^100 kpc of the 
cluster, where most of the excess heating is required. Val- 
ues of / up to 10% might be reasonable if a significant frac- 
tion of the radio plasma seen in cluster cavities with R ^ X 
jBirzan et alJ2004t lG lovannmi & Feretti'206?) is shredded to 
smaller scales, or if the radio plasma producing the diffuse 
cluster radio halo is filamentary rather than homogeneous. 

Thus, even in the absence of microscopic viscosity, a bub- 
bly ICM can prevent leakage of wave energy out of the central 
region. Moreover, bubbles will also boost the dissipation of 
sound waves created by other mechanisms (e.g., cluster col- 
lisions) in the central region, where / is presumably largest. 
Generally, filling factors of / ~ few% — few x 10% are nec- 
essary to provide the required heating in cooling flow clusters. 

4.3. Strong shocks in the cluster center 

Recent deep cluster observations have revealed large 
scale surface brightness discontinuities that are consis- 
tent with modera te strength shocks (McNamara et al. 2005; 
iNulsen et alJl2()05l) . This could mean that cluster centers are 
subject to much stronger shocks than the cool X-ray rims 
around radio bubbles would imply. The question of why 
strong, jet driven shocks have not been found has been a 
puzzle ever since Chandra began observing the temperature 



structures of the inner regions of galaxy clusters. The pres- 
ence of bubbles in the ICM might help explain the lack of ob- 
servational evidence for such shocks: the large sound speed 
inside the bubbles means that they rapidly broadcast the ar- 
rival of a shock, which locally jumps ahead when passing 
over a filament. This not only scatters some of the shock into 
isotropic sound waves, it also distorts and broadens the front 
(Fig. 4). When observing such a shock in X-rays it will appear 
weaker than a coherent shock in the absence of bubbles. 

Such strong shocks might be efficient enough to heat the 
cluster directly. However, the presence of bubbles increases 
the dissipation efficiency regardless of shock strength. There- 
fore, the proposed mechanism would still operate. 

5. CONCLUSIONS 

The presence of bubbles of relativistic gas in galaxy clusters 
can provide a significant boost to the efficiency with which 
shock and sound waves energize the gas. Given the ubiquity 
of both sound waves and of relativistic gas in galaxy clus- 
ters, this may well be a significant (if not dominant) channel 
of cluster gas heating. More generally, our results demon- 
strate the high efficiency of the RMI, which is known to have 
a broad range of scientific and industrial applicability. 
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